The Origin of a Repose Angle: 
Kinetics of Rearrangement for Granular Materials 
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A micro-structural theory of dense granular materials is presented, based on two main ideas. 
Firstly, that macroscopic shear results from activated local rearrangements at a mesoscopic scale. 
Secondly, that the update frequency of microscopic processes is determined by granular tempera- 
ture. In a shear cell, the resulting constitutive equations account for Bagnold's scaling and for the 
existence of a Coulomb criterion of yield. In the case of a granular flow down an inclined plane, 
they accounts for the rheology observed in recent experiments 111 and for temperature and velocity 
profiles measured numerically. |2|,H| Finally, it is shown that for an angle 9 of the plane, the sys- 
tem jams below some critical height H s t op (6), and the critical curve obtained fits remarkably well 
experimental data. 
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An upsurge of interest for granular materials has re- 
cently stirred the physical literature. j^H The problem is 
indeed challenging. The situation is paradoxical. Those 
materials lie at our doorstep, are ubiquitous in every- 
day' life, their understanding is of an extreme interest for 
numberless practical reasons ranging from earthquakes 
or landslides to industrial processes. Yet, there is no 
satisfactory explanation for the obvious fact that heaps 
have a slope. . . Recent advances have been made in ex- 
perimental and numerical studies. In particular, the flow 
of a granular layer down an inclined plane is a laboratory 
model for many realistic situations. ||[i],[7[H ^ n * ms se *~ 
up, evidence has been given for the existence of a critical 
curve H stop (8) relating the angle of the slope, 8, to the 
thickness of a flowing layer below which the system jams. 
This relation refines the well-known Coulomb criterion of 
yield. It is accompanied by Bagnold's scaling, in the bulk 
of a dense flow, which relates the shear stress a to the 
strain-rate e by a oc e 2 . |||| Those experimental findings 
remain unexplained and add even more constraints to the 
challenge faced by the theorist. Kinetic theory [p|JlO|] ac- 
counts for the rheology of dilute systems [jllj , but fails to 
explain jamming and the rheology of dense systems. 

This work focuses on structural rearrangements. It 
draws on the so-called shear transformation zone (STZ) 
theory, [|l2],[13| recently introduced to account for the be- 
havior of viscoplastic solids. The profound difference be- 
tween those two systems, at the microscopic level, is of 
particular interest: it has been pointed out recently that 
some unifying concept might be at work, in various sys- 
tems: granular materials, glasses, foams,... |14| Hence, 
more than an isolated model of granular flows, this work 
attempts to bridge the gap between two of those appar- 
ently different systems, thus showing that structural rear- 
rangement is a key to our understanding of the rheology 
of dense systems on very general grounds. 



Firstly, I will present a basic mechanism for rearrange- 
ment, extracted from [12jlj|. Then, after reviewing the 
basic scaling properties of the TV-body problem, I will 
show how this basic mechanism adapts to the case of 
granular materials. 

A shear transformation zone (STZ) is defined as a lo- 
cus within the material where a rearrangement is made 
possible by the local configuration of the contact net- 
work. [ p2p^ | An important remark that rests at the root 
of STZ theory is that once some microscopic shear has 
occurred somewhere in the material, the system cannot 
shear further at this point, and in this direction (although 
it may shear backward). This leads to the identification 
of pairs of types of arrangements that are transformed 
into one another by a local shear. A local "symmetry" is 
induced by shearing, and the local state of the system is 
determined by the population of arrangements suscepti- 
ble to shear in a given direction. To simplify the picture, 
a single pair of orientations is considered, aligned along 
the principal axes of the stress tensor. Such states are 
schematically represented below: 



although an actual STZ might involve more than four 
grains (say of order 10). 

The populations of STZ's are denoted n + and n_ . And 
R + (rcsp. i?_) denotes the rate at which a + — > — 

(resp. ► +) transformation occurs. Those transitions 

are activated by the free-volume and force fluctuations. 
Since I do not consider variations of volume in this work, 
R± will be written as functions of the stress tensor only. 
Expressions for those rates will be given further. The 
macroscopic shear is seen as the sum of all local rear- 
rangements, 



e = Aq (R+n+ — R-n) 



(1) 



where e is the off-diagonal component of the rate-of- 
deformation tensor, and where Aq is some constant. 

Equations of motion for the populations n± are written 
to close the system. They are of the form, |1^ , ^3[ 

(2) 



h± = R T n T — R±n± + lo (A c — A a n±) 



The first two terms on the rhs account for the "inter- 
nal" dynamics of STZ's, while the last term introduces a 
coupling of the local arrangements with the mean flow. 
This term can be understood in the following way: From 
a macroscopic standpoint, the flow constantly stirs the 
grains, thus creating and destroying local configurations. 
Variable u> is the rate at which the flow induces new con- 
figurations. In a mean-field situation, it can be evaluated 
as the overall work ae produced by the shear stress a, and 
normalized by some typical force of the contact network. 
I will argue in the following that for granular materials, 
this typical force is measured by the pressure P itself, 
so that u; = uj° = aejP if space is not taken into ac- 
count. When spatial extension is considered this term 
introduces a spatial coupling: since an STZ is a config- 
uration of several grains at a mesoscopic scale, the flow 
at a point f contributes to the creation/destruction of 
nearby STZ's. In general, the term u) will be written, 



uj(f) 



1 



w u (r )e ^dr 



(3) 



where a length-scale £ has been introduced, which char- 
acterizes the extension of an STZ. 

The remaining task consists in evaluating the rates R± 
and arguing about the above-mentioned expression for 
u)°. But beforehand, it is necessary to put into perspec- 
tive the essential features and scaling property of granu- 
lar materials. 

It has been observed recently that Bagnold's scaling 
is an immediate consequence of the equations of the TV- 
body problem, so long as no time-scale is imposed to the 
system by coupling it e.g. to a pressure bath. JTl] It is 
however important to understand this scaling in greater 
details so as to write macroscopic equations that agree 
with the underlying physics. 

Consider a system of N grains in contact submitted 
to a set of forces F c at each contact point, and to forces 
F c at contact points with the boundary. The history of 
the system is given by the locations fi(t) of the center 
of mass of each grain, their rotations, the sets c(i) and 
c'(t) of the contact points, and the forces as function of 
time. The interaction between the grains is a pure hard- 
core repulsion; if there is friction, the yield criterion at 
each contact point involves a ratio of forces. Equations of 
motion induce no time-scale, no scale for the forces net- 
work: the system is unchanged if all forces are rescaled 
as, F c ' c — > F c ' c /Fq and if time scales as t —> t\fF~Q. 
Rescaling the forces, leaves the trajectory of the system 
unchanged in the phase space! Only the time-coordinate 
along the trajectory is modified. 



This scaling invariance of the iV-body problem lies at 
the root of Bagnold's scaling. It is noteworthy that ki- 
netic theory correctly incorporates this scaling invariance 
(although it does not account for dense rheology where 
Bagnold's scaling is observed). This is performed via 
the introduction of a granular temperature T which, as 
opposed to thermodynamic temperature, is a dynamical 
quantity, and which determines the collision frequency, 
T = VT. Granular temperature is usually defined as 
the specific kinetic energy. However in a dense mate- 
rial, various contributions can be identified to kinetic en- 
ergy: temperature obviously differs from velocity fluctu- 
ations associated to the rearrangements themselves. In 
this work, granular temperature is seen as associated to 
tiny motion of the grains in the surrounding cage made 
of their neighbors, in the spirit of Haff's fllq , |l6fl picture of 
kinetic theory. By definition, granular temperature de- 
termines the smallest time-scale in the material which is 
also the frequency of all microscopic processes. From a 
phase space perspective, 1/VT determines the time-scale 
at which the system evolves along a trajectory. 

The equation of motion for temperature results from 
an energy balance between the energy produced by the 
flow and energy dissipated by collisions: 



T = ae — aT vT 



(4) 



The second term in the rhs accounts for collision- 
mediated energy dissipation: collisions occur at the fre- 
quency vT, and each collision dissipate an energy pro- 
portional to T itself. Parameter a is related to the resti- 
tution coefficient of the grains. It's straightforward to 
check that this equation in left invariant under the rescal- 
ing a — > a/F , t — ► t^/Fo. This energy balance, added 
to the fact that granular temperature controls the small- 
est time-scale, i.e. the update frequency of microscopic 
processes, is a translation, at the macroscopic level, of 
the scaling invariance of the ./V-body problem. In dense 
flow, the fact that temperature determines the update 
frequency of activated processes, leads to the emergence 
of Bagnold's scaling. 

It is important to make the difference between granu- 
lar temperature and thermodynamic temperature: gran- 
ular temperature does not measure a volume in the phase 
space, henceforth a Boltzmann weight for activated pro- 
cesses. The fact that thermodynamic temperature mea- 
sures a volume in the phase space does not results imme- 
diately from its definition as kinetic energy, by from the 
coupling to a thermal bath, which blends nearby trajecto- 
ries in the phase space, and allows to resort to ensembles. 

Let me now come to the evaluation of the transition 
rates R±. From the preceding discussion, J- determines 
the update frequency. At each trial, the probability that 
grains rearrange depends on the distortion of the force 



network, measured by the ratio a/P. The rates are thus 
written, 

R± = RoVfe ±fMr/p , 

with some constant Rq, and where 1/p is a measure of 
how distorted the force network must be to induce shear- 
ing. There are other ways of understanding why those 
rates should be functions of the ratio a/P: due to the 
hard-sphere interaction, no typical forces is determined 
by some interaction potential, and P appears as the only 
scale of force network; it has also been shown that the 
distribution of forces is of the form exp(— a^Fij < F >), 
indicating that the average normal force is also a measure 
of force fluctuations. |l^ , |l8| 

For similar reasons, the local creation/destruction of 
arrangements performed by the flow, is measured by 
lu° = ea J ' P '. Indeed, creation and destruction rates are 
expected to be proportional to the strain rate. At a given 
deformation, e, the modification of the contact network 
might depend on a/P, but not on any absolute value of 
the forces. Finally, the expression given to ui° can be 
understood as a ratio of the work ai produced by the 
mean flow over some infinitesimal work Pdlo necessary 
to break a contact. 

Let me now write the equations of motion (nh and (0) 
in a more suitable way. Variables 



A 



A = 



n n 



and, ( = 



P 



are introduced, along with the rescaled parameters noo = 
2A C /A a , eo = AoA c /Aa, 7 = AA, and E = 2e R . 
From (|l|) and (||), it comes, 

i = EqVT (A sink (p()- A cosh (fi()) (5) 

A = — (e- 7 wA) (6) 

A = 7 -(l-A) (7) 

with T and u> given by (0) and (0) respectively. Vari- 
ables A and A represent respectively the total normalized 
density of STZ's and the bias between populations n±. 
These state variables account for an history-dependent 
texture (or fabric |0,18[) of the material, and are one of 
the most interesting aspect of the current model. How- 
ever, since only steady state motions are considered in 
this work, A can be safely taken to its asymptotic value, 
A = 1. The fact that A saturates to a constant value 
means that the system adapts to the principal directions 
given by the stress tensor. For the current problems, all 
the complexity of the dynamics lies in the values taken 
by A. 

Before attacking the case of a flowing layer of grains, 
let me consider a simple shear cell experiment, where the 



system admits multiple jammed solutions, e = T = 0, 
for any value of A, as well as a flowing solution for 
A = 1/(7C). In the flowing regime, the dynamics of rear- 
rangements as defined by equations (||) and (g) is stable 
so long as t&nh(fj,() > A. This condition is equivalent to 
C > tan$, where <E> is the solution of, 



tanh(/i tan $) 



1 



7tan$ 



It is a Coulomb criterion of yield with limit angle $. 
Moreover, in the flowing regime, the stationary value of 
the temperature verifies, 



T = 



aE K(() 



with, K(() = sinh(^£) — cosh(p,()/(j<^), whence, 



n3/2 






(8) 



(9) 



which accounts for Bagnold's scaling. 

Let me now consider a flow of granular material down 
an inclined plane. Axis x is taken along the descent and 
axis y perpendicular to the plane. Uniform solutions in 
direction x are looked for, so that the variables e, A, 
and T only depend on the height y. The velocity field 
is oriented along the direction x and the velocity profile 
u(y) is related to the strain rate by, u'(y) ~ 2e(y). The 
total height of the layer of granular material is denoted 
H. The plane makes an angle 6 above the horizontal, 
and the stress tensor obeys Cauchy equations that en- 
force force balance, leading to, P(y) — pg (H — y) cos 9 
and a(y) = pg (H — y) sin#, where p is the mass density 
of the material, and g the gravity. The field ( — tan 9 is 
uniform. 

Let me start with neglected spatial effects, oj = oj°, 
which corresponds to the limit H >> £. The local equa- 
tions are then identical to the mean-field case presented 
earlier. For any y, the systems jams if 9 < $, and flows 
otherwise. In the flowing regime, the temperature and 
strain rate verify equations (|g) and (|9|) respectively, with 
a being a linear function of y. Temperature T grows 
linearly with the depth. After integration over y, equa- 
tion (g) leads to the following expression for the velocity 
field, 



u{y) 



(E K { C)f*^P^(H^-(H-y)y* 
\/a \ 



system is spatially uniform, so that uj 



uJ 



(e. The 



The linear decrease of T with y, this velocity profile as 
well as the rheology u(H) ex iJ 3//2 agree remarkably well 
with experimental and numerical findings. p],p| 

In order to study the effect of the spatial coupling in- 
troduced via uj, it becomes necessary to resort to nu- 
merical integration. The value of A in a steady state is 
now, A = e/'yuj, and varies in space. Typical profiles are 



displayed figure El-(a). The main effects of spatial interac- 
tion appear at the boundaries, resulting in an increase of 
A close to the bottom, while A vanishes at the very top 
(after a small bump due to the convexity of \JH — y near 
y = H) . The increase of A close to the bottom is due to 
the fact that the creation/destruction of STZ result only 
from the action of the flow in the upper half plane y > 0. 
At the top of the flow, A vanishes because e — » 0, whereas 
u> remains non-zero: the system is more liquid-like, while 
it's more textured, solid-like at the bottom. 
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FIG. 1. (a) The state variable A as a function of y/H , 
for values of the parameters taken from (d), angle 8 — 24deg 
and several values of H : H = 5 straight line, H = 10 dashed, 
H — 20 dot-dashed, H — 40 dotted, (b-c) same for the tem- 
perature T, and the velocity profile u. (d) H Btop from experi- 
mental data by Pouliquen fll (circles) compared to the phase 
diagram obtained from the current theory (straight line) with, 
7 = 10.23, fi = 0.98, £ = 0.99. This curve does not depends 
on parameters Eq and a. 

The non-uniformity of A has dramatic consequences 
on the stability of the flow. The jamming criterion 
tanh(y«£) > A is controlled by the maximum value of A, 
reached at the bottom of the flow. Moreover, A globally 
increases with t/H: if the lowest layers jam, the remain- 
ing flow of the upper layers presents an even smaller ef- 
fective height, the criterion for jamming is verified again, 
and jamming propagates upward, resulting in a complete 
arrest of the flow. This process leads to the existence of 
a critical height H s t op {0). Such a curve is displayed fig- 
ure H and compared with experimental data. The fit is 



remarkable. 

The results of this work are two-fold. Firstly, the role 
of activated rearrangements in granular matter has been 
evidenced, hence providing a micro structural interpreta- 
tion for the Coulomb criterion. Secondly, the importance 
of granular temperature as opposed to thermodynamics 
temperature has been shown; the fact that microscopic 
collisions determine the time-scale of activated rearrange- 
ments lies at the root of Bagnold's scaling. Those fea- 
tures contribute to the existence of a jamming criterion 
H s top(0) and to temperature and velocity profiles consis- 
tent with recent numerical results. |3(] 

In any case, the remarkable agreement reached with 
various unexplained features of granular flows, opens up 
new directions in our understanding of dense materials. 
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